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The spectral problem for an integrable system of particles satisfying the fusion rules of 
SU(3)k is expressed in terms of exact inversion identities satisfied by the commuting transfer 
matrices of the integrable fused A^ interaction round a face (IRF) model of Jimbo, Miwa 
and Okado. The identities are proven using local properties of the Boltzmann weights, in 
particular the Yang-Baxter equation and unitarity. They are closely related to the consis¬ 
tency conditions for the construction of eigenvalues obtained in the Separation of Variables 
approach to integrable vertex models. 


I. INTRODUCTION 

Studies of integrable models m in low dimensions have provided important insights into 
the exotic properties of the quasiparticle excitations in a correlated many-body system subject 
to strong quantum fluctuations. Particularly exotic objects are the non-Abelian anyons with un¬ 
conventional quantum statistics expected to be realized in certain fractional quantum Hall states 
GEH: interchanging these quasi-particles can be described by a unitary rotation on the manifold 
of robust degenerate states supported by a collection of a few of them. These degeneracies will be 
lifted by interactions between the anyons and the formation of collective states has been studied 
in models of interacting anyons with given fusion and braiding properties [THEU- 

In this context, the local lattice Hamiltonians generated by the commuting transfer matrix 
of integrable restricted solid on solid (RSOS) models [13J and their generalizations [mfTTj have 
recently attracted renewed interest, see e.g. 0 n ng. In these lattice models the local state 
variables (spins) take values from a given set of representations of a Lie algebra g. The possible 
pairs of spins on adjacent sites are constrained by the fusion rules of the algebra mmm, 
leading e.g. to the RSOS condition for 577(2)*, anyons. Below we consider such models relevant to 
anyons satisfying higher rank fusion rules, specifically of SU(n)y for n = 3. This fusion algebra is 
a truncation of the category of irreducible representations (irreps) of the quantum group U q [su(n)] 
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with q = exp(27ri/(n + k)) or, equivalently, the level-A: integrable representations of the Kac- 
Moody algebra = SU(n ). Local Boltzmann weights for IRF models based on this algebra 

for anyons corresponding to the fundamental vector representation which satisfy a Yang-Baxter 
equation have been constructed by Jimbo, Miwa and Okado m- Using the fusion procedure the 
model has been further generalized to allow for arbitrary SU(n)k anyons relating the admissible 
spins on neighbouring sites in the horizontal and vertical direction, respectively [2D|. Based on the 
fusion hierarchy of functional equations satisfied by the corresponding transfer matrices the central 
charge and conformal weights of the conformal field theories describing the low energy collective 
excitations of the model have been computed mm- 

The goal of this paper is to provide a basis for a different approach towards the solution of the 
spectral problem of this model based on a set exact inversion identities satisfied by the transfer 
matrices (or their eigenvalues) of inhomogeneous generalizations of these models. For the six- 
vertex model and the related RSOS models with integrable periodic or open boundary conditions 
such identities have recently been derived from the underlying Yang-Baxter or reflection equations 
for the local vertex weights using certain physical assumptions such as crossing symmetry and 
unitarity E3H2U, for the vertex models they also arise in Sklyanin’s separation of variables (SoV) 
approach |25]. Complemented with information on the analytical properties of the transfer matrix 
they can be related to the formulation of the spectral problem in the form of Baxter’s TQ-equation 
[26j or inhomogeneous generalizations thereof which arise in models with non-diagonal boundary 
conditions breaking the U( 1) bulk-symmetry of the vertex model [23l [271 |28]. For the six-vertex 
model, the equivalence of these formulations can be shown using the SoV approach where, in 
addition, the completeness of the spectrum has been proven |29j . 

Just as in these cases the Boltzmann weights of the SU(n)k models considered below satisfy a 
unitarity condition. There is, however, no crossing symmetry for the face weights which prevents 
the straightforward extension of the results of Refs. |23l 21] to higher rank symmetries. Motivated 
by relations obtained from the fusion hierarchy of the integrable 5U(3)-invariant vertex model 
and the corresponding SoV [3111 - 132] . we find that a closed set of discrete inversion identities can 
be formulated for two transfer matrices from the fusion hierachy of the IRF model corresponding 
to adjacency conditions in the vertical direction induced by the anyon in the fundamental vector 
representation and its dual, respectively^ 

Our paper is organized as follows: to introduce our notation we first briefly recall the definition of 
the SU(n)k anyon (or A^_ IRF) model and its algebraic structure. Then, motivated by the results 


1 Related identities for SU ( n ) vertex models with various boundary conditions have been constructed in Ref. [33]. 
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for the SU (3) vertex model [30, 3TJ (3U summarized in the appendix, we introduce a generalized 
model based on fused Boltzmann weights and define the transfer matrices appearing in the inversion 
identities for the SU(3)k model. The main result of this paper is the proof of these identities in 


Section IV A The paper ends with a brief discussion. 


II. SU{n) k ANYON MODELS 


Anyonic models can be decomposed into a finite set of topological sectors. The corresponding 
conserved charges obey a commutative and associative fusion algebra 

'tpa <8> i>b - 0 N abi’c (2.1) 

C 

with non-negative integers IVA. In a graphical representaton (to be read from top left to right) of 
fusion the vertex 

b 

a —- c 


may occur provided that IVA ^ 0. The fusion algebra places constraints on the allowed sequence 
of charges in the basis of fusion path states for a model of anyons vp x 


Ipx *Px Ipx 


■ ®n—l®n®n+l ' ) — •• • 


i.e. ipa n+ 1 has to appear in the fusion tp an <g) tp x 


Q"n—1 &n+l 


( 2 . 2 ) 


A. Local states and admissible pairs in the IRF mdel 

For the SU(n)k anyons considered in this paper the sectors are identified with certain irreducible 
representations of the quantum group U q [su(n)} and the fusion algebra is given by the decomposi¬ 
tion of their tensor products into irreps. Just as the level k dominant integral weights of Aj l _ 1 the 
topological sectors in the SU(n)k anyon models are represented by vectors 

n —1 n —1 

a = ai = k (2-3) 

2—0 2—0 

with nonnegative integers ai and being the fundamental weights of A^ 1 2 1 with uj n = uq. For 

given k > 1 the topological charges (or the spin variables in the IRF model) take values in the set 
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[0,0] 

[1,0] 

[hi] 

[2,0] 

[2,1] 

[2,2] 

[1.0] 

[1,0] 

[1,1] © [2,0] 

[0,0]® [2,1] 

[2,1] 

[1,0] ©[2, 2] 

[1,1] 

[1,1] 

[1,1] 

[0,0] ©[2,1] 

[1,0]® [2,2] 

[ 1 , 0 ] 

[1,1] © [2,0] 

[2,1] 


TABLE I. Fusion rules for SU ( 3)2 anyons involving the fundamental anyons and their adjoints corresponding 
to the Young diagrams [1,0] and [1,1], respectively. 


P+(n; k) of dominant weights (2.3). A convenient representation of these local states is in terms of 
Young diagrams [Ai,..., A n _i] with A i nodes in the z th row: for k = Aq > Ai > ... > A n _i > \ n = 0 


(2.4) 


this diagram can be identified with the local state 

n —1 n —1 

g = (At - Aj+i) Ui = A0W0 + ^ A* e, G P+(n; k ) 

2—0 2=1 

with the elementary vectors e* = w* — 1 . Two diagrams represent the same element if and only 

if one is obtained from the other by removal of columns of height n. 


As discussed above the fusion algebra (2.1) leads to a set of constraints for an ordered pair of 


local states (a, b) with a, b e P + (n; k) to be admissible as a configuration on neighbouring sites of 
the lattice model. For a system of anyons ij)yx\ corresponding to the Young diagram [A] such a pair is 
allowed if ^ 0. Considering anyons in the fundamental vector representation [1] = [1, 0,..., 0] 
the vertex 


may occur in a fusion path state provided that 


b = a + e i , for some i = 1,..., n , 


(2.5) 


see e.g. Table [T] for SU( 3)2 anyons. 

The set of local states and the constraints for a given type of anyons can be represented in 
an oriented graph with nodes labeled by elements of -P+(n; k) and arrows from node a to b repre¬ 
senting an admissible pair (a, b). The underlying fusion rules for 'i/im anyons are encoded in the 
corresponding adjacency matrix with elements 


ab 



( 2 . 6 ) 


Note that for [A] being a symmetric tensor or antisymmetric tensor the decomposition ?/> a <8> ip[x] is 
multiplicity free and the elements of the adjacency matrix are 0 or 1. The adjacency graph for the 
local states P + (3; 2) in a system of 7/ip] anyons is shown in Figure [lj 
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FIG. 1. Adjacency graph for the set of local states P + (3; 2). Admissible pairs of states in the anyon 
model corresponding to the fundamental vector representation [1] are denoted by an arrow connecting the 
corresponding vertices. For anyons f/’fi 2 ] hr the adjoint representation the arrows would have to be reversed. 

B. Boltzmann weights 

The j IRF model for SU(ii)k anyons in the vector representation [1] is defined on a square 

lattice, such that the spins on the corners of a face take values in P+(n;/c). A Boltzmann weight 

corresponding to a configuration a,b,c,d round a face is depicted as 

x c d 

c d 

' u I = 


W 


b a 


u 


b a 

It is non-vanishing if and only if the four pairs of variables (a, b), (6, c), (d, c), (a, d) are all admissible 


in the sense of (2.5). Their explicit expressions read |T5 

a + 2e* a + ei 

CL Ci CL 

a + e* + ej a + ej 


W 


IT 


W 


a + ei 


CL €i Cj CL Cj 


a + ei 


u ) = n (u) [u + 1] 

U + a ji\[ 1 ] 


u = n(u) 


u = n (u) 


[' a ji\ 

M[°ji ~ !] 

[ a ji\ 


for i ^ j 


for i ^ j . 


(2.7) 


For a state a 6 -P+(n; k) given by Eq. (2.4), is defined as the inner product 

aij = {a + p+i- ej) = i + (2.8) 

with p being the sum of fundamental weights. At criticality, the dependence of the Boltzmann 
weights on the spectral parameter u is given by trigonometric functions 

7r 


[u] = sin(7i7y), ij = 


n + k 


(2.9) 
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These weights satisfy the initial condition 



( 2 . 10 ) 


Notice also that at the special value u = 
a —> a + 2ej, are automatically excluded, 
equation (YBE) 


— 1 the weights containing straight paths, that is from 


The Boltzmann weights (2.7) satisfy the Yang-Baxter 


Vw 


a g 

f e 


Vw 


v W 


a b 
f 9 


a b 


9 c 


u W 


g c 

e d 


u + v 


u + v\W 


f g 

e d 


u W 


b c 

g d 


( 2 . 11 ) 


and an inversion relation which can be used to fix the normalization function r\(u) 


E w 


d g 
a b 


u W 


d c 

g b 


-u\ = d ac rl(u)[l - u][l + u\ 


( 2 . 12 ) 


There are additional symmetries, e.g. symmetrizability, Z n -invariance and a duality, see [15] for 
details. 

For the proof of the inversion identities below we use graphical representation of these local 
relations. In particular, the YBE may be depicted as 


b b 



u 


U + V 


C / 


'g 

U + V 

9 

u 



d d 


In this and in the following figures, nodes marked with a dot (•) represent spin variables which are 
summed over all possible local states and nodes with equal spins are connected with a dotted line. 
Unitarity condition is depicted in a similar fashion as 
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C. Hecke algebra and projectors 


The algebraic structure underlying the A^_ , RSOS models is connected with a quotient of the 
braid group, namely the Hecke algebra [2D ESI ESI- This connection leads to the construction of 
projection operators, which are used later in the text in order to prove the inversion identities 
satisfied by the transfer matrices. 


Let first |ao • • • ol+i) be a fusion path state (2.2) for the anyons considerered, i.e. a sequence of 


local states a* such that each pair (oi,aj + i) is admissible. We define the Yang-Baxter operators by 
their action on these states as 


W i (u)\a 0 ---a L+1 )= W 

bieP+(n-k) 


( tti—1 

.) 

y b{ o>i +1 

) 


U I \dQ • • • Cli— \biCLi -\-1 • • • 


(2.13) 


for i = 1,..., L. In terms of the Yang-Baxter operators, the YBE (2.11) takes the form 


Wi(u) W i+l {u + v) Wi{v) = W i+ 1 (v) Wi(u + v) W i+ 1 {u ), 


(2.14) 


which is reminiscent of the braid relation satisfied by the generators of the braid group. In fact, 


by proper choice of the normalization r\(u) of the Boltzmann weights (2.7) these operators can be 
written as 


u 


Wi(u) = e~ ir)U I + e~ l1) Xi. 

Here the Xi s are independent of the spectral parameter and satisfy 


(2.15) 


Xi X i+ i Xi — Xi + 1 Xi Xi + 1 

X? -{q-l)Xi~q = 0 (2.16) 

[X i ,X j \=0, \i — j\ > 1 

with q = exp(2 ii]) = exp(27r i/(n + k)) being the deformation parameter of the underlying quantum 
group U q [su{n )]. This set of relations gives rise to a representation of the Hecke algebra HL+i(q)- 
Note that in the n = 2 case, that is in the RSOS models, the Hecke algebra truncates essentially 
to the Temperley-Lieb algebra. As was pointed out in [2D one actually obtains the representation 
of the quotient of the Hecke algebra in which g-analogues of the full Young (anti-)symmetrizers 

(P^ 1 ) 2 = Pf for the SU{n)k model can be defined. They act on (admissible) sequences of local 

states of length i + 2 and satisfy 


P n (*, ■ ■ •, * + n + 1) = 0 . 
(*,...,* + 1 + 1) = 0 , 


Vi, £ = k + l, ...,n + k— 1. 


(2.17) 
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In terms of the Yang-Baxter operators they can be formally written as 


l -1 

m =0 



[i] 

[j + 1 ] 


Wj(±j) 


(2.18) 


As an example, we have the following schematic representation for P 1 2 (as before the spin variables 
Cj on dotted nodes are being summed over) 



The action of the projectors can be deduced from the explicit expressions of the Boltzmann weights 


(2.7). First, one may observe that for u = — 1, the first of the weights (2.7) is always zero. Hence the 
operator Pj - projects out all states ( 02 , do) which are connected by a straight line in the adjacency 
graph of local states P+(n; k), see Figure [IJ Similarly, it can be shown that P~_i|ao ■ ■ • a n ) oc <5 an a 0 
for any fusion path state |ao • • • a n ) of SU(n)k anyons. 


III. FUSED A^2 1 IRF MODELS 
A. Commuting transfer matrices 


Starting with the set of Boltzmann weights (2.7) satisfying the Yang-Baxter equation we can 
define a family of transfer matrices generating the commuting integrals for a system of SU(n)k 
anyons in the vector representation. Introducing local inhomogeneities {ufc}^ =1 the first transfer 
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matrix for the A^}_ 1 IRF model is constructed as usual 


t(u)= n w 

&k—1 &k 

u - u k 

k=1 

y bk-1 bk 

) 


a o a i Ofc-i fflfc cll- 1 aL 


| U — Ul 


u — u k 


, u — ul 

(3.1) 

b 0 ' 

>i bk 

>L -1 &L 


after identifying (az,, bi) = (ao, bo) for periodic boundary conditions. The quantum Hamiltonian for 
the lattice model of SU(n)k anyons VR] with local interactions is obtained from the homogeneous 
limit Uk —> 0 for all k of the transfer matrix and can be expressed in terms of the Yang-Baxter 


operators (2.13) 


H = cUnT(u)| u=0 = d u \nWi(u)\ 


u =0 


(3.2) 


2=1 


To obtain the complete set of commuting integrals one has to consider fused A^}_ ] IRF models, 
similar as in the case of the integrable SU(ri) vertex models (30 ) 134 j . see Appendix |A| for n = 3. 


Analogous to (A7) we define the second transfer matrix from this family of operators acting on the 


space of fusion path states of SU(n)k anyons in the horizontal direction by 


U(u) = Pf T(u) T(u - 1) 
ao bo 
bj \ bL Co 




1 n w 


k=1 


O'k—l &k 
bk-i bk 


u-Uk \ W 



U — Uk — 1 


(3.3) 


a i 


dk- 1 


dk 


a-L -1 


a L 



12 — 221 


U — Uk 




U — Mfc — 1 



u — u L 


b L 


u — Ul — 1 


Cl 


Ck-1 


Ck 


cl -i 


cl 


where, again, (cil,cl) = (aoTo) for periodic boundary conditions. Further transfer matrices are 
written as 


T e (u) = P^T^Tiu - 1)... T(u - £ + 1), 


(3.4) 


for £ = 3,... ,n. Note that P n _ 1 projects on a one-dimensional space, therefore T n (u) = A(w)I, 


similar as the quantum determinant (A8) of the vertex model. The function A {u) depends on the 
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type of SU(n)k anyons considered, for the present case of we find 


A (u) 


L / n— 1 

n I [u - u k + 1 ] Y\[u-u k -1] 
k= 1 V 1=1 


(3.5) 


B. Fusion of weights 

Alternatively, the fused , IRF models can be constructed by means of the fusion procedure 
. Boltzmann weights VF^^u) with admissible pairs corresponding to fusion with the anyon 
V’fA] (^[/i]) along horizontal (vertical) links are obtained from the partition function for a suitably 


chosen rectangular block of the elementary ones (2.7) with properly shifted spectral parameters 


by projection onto the Young diagram [A] ([/r]). In the present context where the admissible 
fusion path states in the horizontal direction are those of the SU(n)k anyons '(/'[I] we can restrict 
ourselves to the vertically fused Boltzmann weights W^^ L \u). Specifically, we consider the case of 
[/i] = [l 2 ] = [1,1,0 ,...] which can be obtained by fusion of two elementary weights. 

Following the prescription used in Ref. m the fused Boltzmann weights are constructed graph¬ 
ically as follows (fused edges are represented by double lines on the corresponding link, as before 
arrows indicate the direction in which the constraint for admissible pairs is to be read) 


w wi 1 


2i / 02 i>2 


ao bo 


u \ =W V 


02 


a o 


a 2 b 2 
oo bo 


a 2 



b 2 




-'ll 


Here one has to perform an antisymmetrization on the nodes marked by ©, in addition to summa¬ 
tion over all possible values, corresponding to the projection onto [l 2 ]. The fused weights do not 
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depend on the spin b\. They can be written in the compact form 

2 


W v 


W v 


w v 


a T ej, + e\ a + e\ 


a + e 


I’k 


a 


u = - 


r 2 {u)\u- £][« + 5 ] U 
3 = 1 






a. + + eA a + eA 

a + e im a 

a + e a + ej m a + eA 
O + a 


U = -r 2 (u)[u- \][u+\ +a ikim ] 


[ 1 ] 


n 


- 1 ] 


' l k' l m J j = \ 

j^m 




(3.6) 


u = -r 2 (u)[u - \}{u+ |] J] 


l=i 


[ a b*m 1] 


where we have defined e\ = eq + ej 2 and require that m G {1,2} and ik ^ {* 1 , * 2 }- The spectral 
parameter dependent scalar factor r 2 {u ) is fixed by our choice of normalization of the elementary 


weights (2.7), see m- 

The antisynmretrization of the vertically fused weights allows to write them as a column of 


two weights, multiplied by the projector P 1 , Eq. (2.18), excluding straight paths in the adjacency 
graph. In terms of a graphical representation, we have 


02 


a 0 



b 2 

h 

b 0 


This equivalence will turn out to be particularly useful when we consider the transfer matrices and 
their algebraic relations in the next section. 

As a final note, it should be stressed out that these fused weights satisfy a set of Yang-Baxter 
relations, ensuring the preservation of integrability. In Figure [2] we represent the Yang-Baxter 
relations satisfied by the elementary and fused Boltzmann weights in addition to the initial one 


( 2 . 11 ). 


In terms of the vertically fused Boltzmann weights the transfer matrix (3.3) is written as 


do 


a 1 


Ofe-l 




a,L -1 


dL 




1 


1 

|W — Ml — 2 


U Uk 2 

‘ ‘ ‘ 

U — ILL — 2 ‘ 



Cfc-1 


Ck 


cl- 1 


cl 


(3.7) 


11 Uk 


k= 1 


relating it to the fused } IRF model based on fusion paths of and ^[ 12 ] anyons in the 
horizontal and vertical direction, respectively. Similarly, the transfer matrices Te(u) are obtained 
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FIG. 2. The additional Yang-Baxter relations containing the fused Boltzmann weights (3.61. Double lines 


represent fusion along the particular direction. Indices and spectral parameters are suppressed for the sake 
of clarity. 

from the IRF model with ^[l^] anyons on the vertical links. We do not consider transfer matrices 
with anyons other than those related to these fundamental representations on the vertical links here. 
Based on our results below we conjecture that they are algebraic functions of the T e(u), similar as 
in the integrable SU{n ) vertex models [5Uj f33], thus do not provide additional information on the 
system. 


Note that as a consequence of (3.6) we have U(itfc + 1) = 0 for k = 1 ,L, resembling Eq. (A12) 
for the vertex model. This allows to extract a spectral parameter dependent factor from the second 
transfer matrix 


u (u) = ( IJ [« - u k ~ 1] ) u (u). 


(3.8) 


^fc=i 


IV. INVERSION IDENTITIES FOR THE a! 1} IRF MODEL 


Motivated by the results for the SU( 3) vertex model given in the Appendix, our aim is to 
transfer them to the face model. Initially, the or RSOS models have been introduced by 
Baxter in order to solve the eight-vertex model |38H40| . The observed correspondence between 
Boltzmann weights of faces in the RSOS models and vertices in the eight-vertex model turns out 
to be quite generic and can be formulated based on the identification of equivalent operators in 
face and vertex models, respectively [4T|. For the A^ 1 _ 1 IRF models the face/vertex correspondence 
was stated in Ref. M- 

As we have seen above many of the quantities arising from the integrable structures in the SU (n) 
vertex models can be defined similarly in the (fused) IRF models. Inspired by this observation 


and the identities (A13) satisfied by the transfer matrices of the SU(3) vertex model we now show 


hat similar relations hold for the transfer matrices T(rt) and U(tt) of the inhomogeneous IRF 































models, namely: 


T (u k ) T (u k - 1) = U (u k ), 

T (u k ) U (u k - 1) = A (u k ), 

U(uk) U(n fc - 1) = A (u k ) T (u k - 1) 
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(4.1) 


for k = 1,... ,L. Not all of these identities are independent: as in the vertex model each relation 
for given k can be obtained from the other two. Furthermore, due to the translational invariance 
of the model, one has e.g. 

L L 

n T (“*)= ( n [«fc-«<+i])i- (4.2) 

k=l k,£= 1 

For the proof of these identities we make use of the properties of the Boltzmann weights in¬ 
troduced above adapted to the rank n — 1 = 2 case: First, we note that for the case of the 
model the anyon ^ 2 ] corresponds to the adjoint of the vector representation. Therefore the anti- 
symmetrization in the fusion of Boltzmann weights effectively reverses the order of local states in 
the admissible pairs or, graphically, 


This can be used to show that, upon proper normalization of the W v , we have in particular 

d 


c = $bd Cab [u + 1] [u - 1] [u - 2] , 



where the coefficients C are given by 


Cab — ( 


if a is a corner state in the adjacency graph 
f2l 

Tjj if b is a corner state in the adjacency graph 
1 otherwise. 


(4.3) 


(4.4) 


Furthermore, since Ylk=i e n = 0, the definition (3.6) of the fused weights implies that, for n = 3, 


they satisfy the initial condition for u = — | 

W„ 


fed 

A 

a 

~v 


5an ■ 




(4.5) 
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Note that a similar relation holds for the Boltzmann weights obtained by fusing to — 1 faces in the 
SU(n)k anyon model. 


A. Proof of the inversion identities 


The first of the inversion identities (4.1) follows trivially from the definition (3.3) and the fact 


that the k-th column of T{u k )T{u k — 1) reduces to the projection operator P 1 . Using the Yang- 
Baxter equation (2.11) and {Pi) 2 = P^~ yields the desired result. 


To prove the second of the inversion identities (4.1), it is useful to employ the expression of U 


in terms of the vertically fused weights (3.7) and exploit the inversion relation (4.3) satisfied by 


them. Define also Ujj = m — Uj. We have then 


T(u k )XJ{u k - 1) = 


Oo a 1 flfe-l CLk CLL-l O.L 


Ukl 


0 

ibfe-i , , 


UkL 

1 Ukl - I 


3 

2 


UkL § 1 


Co ci Ck-1 Ck cl-i cl 


Directly from the definitions of the elementary Boltzmann weights (2.7) and the vertically fused 


ones (3.6), it follows that the above expression contains a factor 5 ak ^ k _ 1 5 Ck ^ k _ x ~ ba k c k - Repeated 


use of the inversion relation (4.3) between W v and W then leads to 


T(itfc) 1) | [u k £ T 1] [u k p 1] [u k £ 2] C 

f 


cece +i 


1=1 


(4.6) 


= A(u fc )I, 


where we have used the fact that, as a consequence of periodic boundary conditions for our model, 


JJ Cc t c l+1 — 1 • 


(4.7) 


t =l 


Alternatively, the proof can be done using the projector operators discussed in Section IIC (in fact, 


this form will turn out to be more convenient in order to prove the third of the inversion identities 
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below): 

T(u k )XJ(u k - 1) = 


a o a l fflfc-i a,k cll -i cll 




In the last step we have repetitively used the Yang-Baxter relation in order to push the left projector 
faces, to the right, and exploited once again the fact that (Pf) 2 = Pf. One immediately 
recognizes then the formation of the second projector operator, Pf, constructed in Section IIC 


As already noticed in that section, the action of P 2 ~ on a sequence (ao,... , 0 , 3 ) is proportional 
to 5 ao a 3 , the exact factor of proportionality being the one derived above through the use of the 
vertically fused weights. 
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For the third of Eqs. (4.1) we have in the same spirit 


U(uk)U(u k - 1 ) = 


do di d;c_i ak ciL-i a L 



YBE 



b L 


d-L-i 


YBE 


do di dfc^i 



dk dL -i dL 
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YBE 


ao ai a fe _i 



dk dL -1 < 4 l 


YBE 


ao ai dk-i 



inv. id. 


A(u k ) 


= A(u k ) 


a o ai ak -1 







dk 

dL 

-1 

dL 

ao 

ai 

dk-i 

dk 

dL-1 

dL 


Ukl — 1 


Uk Uk 1 1 

UkL — 1 1 

do 

di 

dk- 1 

dk 

dL- 1 

dL 


= A (uk) T (u k - 1) . 


V. DISCUSSION 


In this paper we have proven a set of discrete inversion identities (4.1) satisfied by transfer 


matrices of inhomogeneous versions of an SU(3)k anyon chain (or IRF model). As shown 
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in the appendix, related identities can be derived using the quantum inverse scattering method 
(QISM) for the SU(3) quantum spin chain (or vertex model). 

Similar identities have been obtained for the six-vertex model within Sklyanin’s separation of 
variables solution and for the RSOS models starting from local properties of the corresponding 
Boltzmann weights [241 [25 j. In these cases, both with underlying rank-1 quantum group, the in¬ 
version identities, when complemented by information on the analytical properties of the transfer 
matrix, have been found to provide a formulation of the spectral problem which allows to compute 
all eigenvalues. We emphasize, however, that such purely functional approach has to be comple¬ 
mented by an independent check that a given solution actually corresponds to an eigenvalue of the 
quantum chain. For the six-vertex model the latter is provided by the SoV approach. 

To what extend similar results can be established for the rank n — 1 = 2 models considered here 
remains to be studied: for the SU( 3) vertex model some evidence exists from the SoV approach 
m but various open questions, e.g. concerning the actual construction of the separated variables 
in the fundamental spin representation (which is necessary to compute eigenstates) and the nature 
of their common spectrum, remain to be addressed. In addition, eigenvalues obtained within the 
functional approach for the vertex model can be checked, for the periodic boundary conditions 
considered here, against those from the algebraic Bethe ansatz m- 

For the IRF model we have some preliminary numerical results for small systems, but there 


remains the practical issue that identities such as (4.1) (or (A13) for the vertex model) are not 
well suited for an efficient computation of the transfer matrix eigenvalues. Reversing the line of 
arguments used to obtain the inversion identities in the appendix, however, they can be related 


to generalized TQ-equations such as (All) for the SU(3) vertex model. This requires to find a 


factorization of (3.5) compatible with the asymptotic behaviour or the transfer matrices. Part of 
the additional input required to address this question for the critical IRF models is available: by 
definition the transfer matrices appearing in the inversion identities are Fourier polynomials in the 
spectral parameter. The underlying fusion algebra allows to split the spectrum into topological 
sectors m where the asymptotics of the transfer matrix can given in terms of the eigenvalues of the 


adjacency matrix (2.6). As a consequence, the effect of the anyonic statistics on the spectrum of the 
IRF model as compared to the vertex case is similar to that of a twist in the boundary conditions 
(depending on the sector), in agreement with the results for the An * IRF models obtained from 
the fusion procedure [21). 

Finally let us note that - while we have concentrated in this paper on the derivation of inver¬ 
sion identities for the the critical rank-2 IRF models with generic inhomogeneities and subject to 
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periodic boundary conditions - we expect that similar identities can be derived for the transfer 
matrices T e(u), t = 1 ,,n — 1, for the A^_i IRF model and for models with open boundary 
conditions }42] - similar as in the case of the SU (n) vertex models {32]. 

We plan to address some of these questions in future work. 
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Appendix A: The SU{ 3) vertex model 

In this appendix we recall the structures underlying the integrability of the SU( 3) invariant 
quantum spin chains and related two-dimensional vertex models and their Bethe ansatz solutions 
[3oi eh m S3]. The Hilbert space of the vertex models is the tensor product T~L = <8>^ =1 Vj. 
To be specific we consider the case where Vj ~ C 3 is the space of quantum states at site j of 
the lattice corresponding to the fundamental three-dimensional (vector) representation of 517(3) 
(corresponding to the Young diagram [1] = [1,0]). In the framework of the quantum inverse 
scattering method (QISM) we define the monodromy matrix acting on the tensor product of the 
auxiliary space V a ~ C 3 and the Hilbert space R of the model as 

Ta(u) = C aL (u - U L )C a ,L-l(,U - U L - 1) ■■■ Cal{u- U\) . (Al) 

Here the C a j ( u ) are operators acting non-trivially only on V a ®Vj. They are given as 

3 

Caj (u) = U I ( 8 ) I + ^2 e u® e £k ( A2 ) 

k,e =i 

with 3x3 matrices (= Skm.Sen acting on V a . The complex parameters Uj. j = 1,..., L, 

V / mn 

define inhomogeneities in the lattice. 

The monodromy matrix (as well as the local /^-operators) is a representation of the Yangian 

y(su( 3)) 

Rab(u - v)Ta{u)Tb{v) = Tb{v)T a {u)R ab {u - V) . (A3) 

The R-matrix comprises the structure constants of this algebra 


Rab(u) = it I (8)1+ P a fc, 


(A4) 
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where P a {, is the permutation operator acting on the tensor product C 3 <8> C 3 as P a j,(x ® y) = 


y ® x , Vx, y E C 3 . As a consequence of the Yang-Baxter relations (A3) the transfer matrix 


T(u) = tr a [T a (u)} 


(A5) 


forms a family of commuting operators, [T(u ), T(v)] = 0. Starting from the reference state 


®j = i |0) ■, where each spin is in the highest or lowest weight state of the local SU(3) irrep, the 
spectrum of the transfer matrices can be obtained by means of the nested (coordinate or) algebraic 
Bethe ansatz (ABA) [Ml 03 33] . 


An alternative solution of the spectral problem relies on functional relations between (A5) and 


more general SU (3) symmetric transfer matrices acting on the Hilbert space T~L of the vertex model 
. With the projectors 


1 


1 


^ab — I^ab) — 1) 

Fabc = + PafePfoc + P 


- Pfec) , 


(A6) 


onto the antisymmetric subspaces of the product V a ® Vj, and V a ® V& ® V c , V a ~ C 3 , respectively, 
we define a second transfer matrix U (it) 


U(u) = tr ab [P ab T a {u - l)7i(u)] 


(A7) 


generated by £-operators corresponding to the adjoint [l 2 ] of the vector representation of SU(3) in 
auxiliary space V a — C 3 and, similarly, the quantum determinant of the monodromy matrix T{u ), 


A (it) = t^abc[P abc Ta{u - 2)T b {u - l)Tc(u)\ 


(AS) 


which generates the center of the Yangian y(su(3)). By construction A(it) is a polynomial in the 
spectral parameter. It has c-number valued coefficients and can be factorized as 


A(it) = di(it — 2 )d,2(u — l)ds(u) I. 


(A9) 


The polynomials dj{u ) depend on the representation of the Yangian in question. Here, i.e. for the 


inhomogeneous model (Al) based on the vector representation of SU{ 3) in all components of the 


quantum space, they are found to be m 


L L 

di(u) = JJ(ii - Uj) = d 2 (u ), d 3 (u) = JJ(it - uj + 1). 

3 =1 3 =1 


(A10) 
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The transfer matrices T(u ) and U[u) generate the complete set of commuting integrals of the 
SU( 3) spin chain. They satisfy functional equations with auxiliary operators Qi :2 (u) 


d 2 {u — 2) d^u — 1) Q\(u — 3) — U(u - l)Q\(u — 2) 

+ d\(u — 2) T(u — l)Qi(u — 1) — d\{u — 1) d\(u — 2) Q\(u) = 0 . 
d 3 (u - 2) d 3 (u - 1) Q 2 (w - 3) - d 3 (u - 1) T(u - 2)Q 2 (u - 2) 

+ U(u — l)Q 2 (u - 1) - di(u - 2) d 2 (u — 1) Q 2 (u) = 0, 


(All) 


similar to Baxter’s TQ-equations for the transfer matrix of the eight-vertex model [26] . As a 
consequence of the commutativity of the transfer matrices and the Q-operators among each other 
for different arguments, analogous third order difference equations holds for their corresponding 
eigenvalues. Using the analytical properties of the transfer matrices these eigenvalues can be 
computed reproducing the result obtained from the ABA [30]. 

Note that the actual solution of the spectral problem for the transfer matrices by means the 
Bethe ansatz methods introduced so far is not possible for all integrable lattice models: the ABA 
relies on the existence of a suitable (highest or lowest weight) reference state which does not always 
exist, e.g. for models with boundary conditions breaking all possible 17(1) symmetries. Similarly, 


the functional approach based on the TQ-equations (All) requires a sufficiently simple (e.g. poly¬ 


nomial) parametrization of the eigenvalues of the Q-operators. Neither of these requirements is 
met, e.g., for spin chains with open boundary conditions subject to non-diagonal boundary fields. 
For models with U q [.s > u(2)]-symmetry this issue has been addressed recently using separation of 
variables (SoV) and through the derivation of inversion identities satisfied by the transfer matrices 
for certain arguments, both leading to certain generalizations of the TQ-equations |24[ [28l [29l 3S] . 
An added value of the formulation of the spectral problem within the SoV approach is that it 
provides a basis for the proof that the solution is complete 123 [29J S 3 HE]. 

Application of Sklyanin’s SoV approach to integrable SU (3) models leads to equations similar 
to the TQ-equations but with Q \ ^-eigenvalues being functions on the discrete set of common 
eigenvalues of the separated coordinates m- By explicit construction for small systems we find 
that, similar as in the SU(2) case, this set is contained in the integer spaced lattice of u -values 


enclosed by the singular points of the difference equations (All). For the model with local spins 


carrying the fundamental representation, i.e. (A10), the spectral parameter u takes values from 
{uk,Uk + 1, Uk + 2}(j =| . Eliminating the corresponding amplitudes Qi, 2 (^) one finds that 


U(u k + 1 ) — 0 . 


(A12) 
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and arrives at the following set of inversion identities for the transfer matrices 

T(u k )T(u k - 1) = U(u k ), 

T(u k ) U(u k - 1) = A (u k ), (A13) 


U(u k ) U(u k - 1) = A (u k ) T{u k - 1), 


for k = 1,..., L. Using the projection property (A6) of the .R-matrix and the Yang-Baxter relations 
satisfied by the transfer matrices T(u) and U (u) these product identities have been derived before 
by Cao et al. [33] (note that the third identity can be obtained from the two other ones). These 
equations, together with the analytical properties of the transfer matrix are sufficient to compute 
their eigenvalues. 

We end this appendix by noting that — while we have considered the Sf7(3)-invariant rational 
vertex model - it is straightforward to extend the discussion to the anisotropic (g-deformed) model 
with trigonometric dependence of the vertex weights on the spectral parameter [35] SMU. 
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